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ABSTRACT 

Substantial orbital migration of massive planets may occur in most extrasolar planetary sys¬ 
tems. Since migration is likely to occur after a significant fraction of the dust has been locked 
up into planetesimals, ubiquitous migration could reduce the probability of forming terrestrial 
planets at radii of the order of 1 au. Using a simple time dependent model for the evolution of 
gas and solids in the disk, I show that replenishment of solid material in the inner disk, following 
the inward passage of a giant planet, is generally inefficient. Unless the timescale for diffusion of 
dust is much shorter than the viscous timescale, or planetesimal formation is surprisingly slow, 
the surface density of planetesimals at 1 au will typically be depleted by one to two orders of 
magnitude following giant planet migration. Conceivably, terrestrial planets may exist only in 
a modest fraction of systems where a single generation of massive planets formed and did not 
migrate significantly. 

Subject headings: accretion, accretion disks — planets and satellites: formation — solar system: 
formation - planetary systems: formation — planetary systems: protoplanetary disks 


1. Introduction 

Radial velocity surveys (Mayor & Queloz 1995; Marcy & Butler 2000) have discovered a significant 
population of giant extrasolar planets with orbits interior to the snow line of protoplanetary disks (Hayashi 
1981; Bell et al. 1997; Sasselov & Lecar 2000). These discoveries - together with the generally accepted 
theoretical belief that giant planet formation is difficult at such small radii (Bodenheimer, Hubickyj & 
Lissauer 2000) - suggest that massive planet formation is often accompanied by orbital migration, which can 
be driven by gravitational interactions with the gas disk (Goldreich & Tremaine 1980; Lin & Papaloizou 1986; 
Lin, Bodenheimer & Richardson 1996). The ss 2/3 of currently known systems that have semi-major axis 
a < 1.5 au are obviously prime candidates for having undergone migration, but theoretical models (Trilling 
et al. 1998) suggest that even this is an underestimate. Trilling, Lunine & Benz (2002) have shown that 
the observed planets were probably preceded by earlier generations that failed to survive migration, while 
Armitage et al. (2002) estimate that little or no migration (defined as a fractional change in semi-major axis 
|Ao|/a < 0.1) occurs in as few as 15% of systems with a surviving massive planet. 

If migration is indeed ubiquitous, it is reasonable to expect that it will affect both the formation and the 
survival prospects of co-existing terrestrial planets. Several authors have studied the dynamical stability of 
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notional terrestrial planets in known extrasolar planetary systems (Gelnnan, Adams & Laughlin 1996; Jones, 
Sleep & Chambers 2001; Menou & Tabachnik 2002). Stable orbits that could in principle harbor terrestrial 
planets exist in most of these systems. The exhaustive study by Menou & Tabachnik (2002) concludes that 
at least 25% could host habitable terrestrial planets, while uninhabitable low mass planets could be present 
in a larger fraction of systems. 

In this paper, I investigate whether migration is likely to restrict the formation of terrestrial planets. 
It is fairly clear that a doomed giant planet, migrating inward toward eventual merger with the star, will 
cause any pre-existing low mass planets or planetesimals at smaller radii to be lost. What is less obvious is 
whether a subsequent generation of planetesimals (and eventually planets) could form from the potentially 
still massive disk present after the merger. To investigate this, I describe in §2 a simple model for the 
evolution of gas and solids in the disk. In §3, I present results for the surface density of solids at small radii 
following massive planet migration. I find that rather little material is likely to be able to repopulate the 
terrestrial planet forming region, and discuss in §4 the implications for the frequency of terrestrial planets. 


2. Gas disk and dust evolution model 


The evolution of the protoplanetary disk, and the coagulation of dust within it to form (eventually) 
planets, remain imperfectly understood. I describe here a simple three component model, for the evolution 
of gas, dust and planetesimals, that is intended only to capture the broad outlines of the probable behavior. 


The gas component is modeled as a viscous accretion disk (Pringle 1981), with disk parameters chosen to 
reproduce the typical lifetimes (Strom et al. 1989; Haisch, Lada & Lada 2001) and accretion rates (Gullbring 
et al. 1998) of Classical T Tauri disks. The surface density E(r, t) evolves according to, 
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where v is the kinematic viscosity and the term E w i n d allows for mass loss in a disk wind. For the viscosity 
I adopt, 

!/ = 4.1 x 10 13 (-!—) cm 2 s" 1 , (2) 

V1 au/ 

while for the mass loss a prescription motivated by observations and theoretical modeling of photoevapora¬ 
tion (Johnstone, Hollenbach & Bally 1998; Clarke, Gendrin & Sotomayor 2001; Matsuyama, Johnstone & 
Hartmann 2002) is used, 


Ewind CX r 5 / 2 , r > 5 au 

E w ind = 0; r < 5 au. (3) 

The normalization of the mass loss is chosen such that the mass loss rate within 25 au is 10~ 9 M 0 yr , 
The rate of surface density loss E w ind is constant in time. The initial surface density profile is, 



where ri n = 0.067 au is the inner disk radius and ro = 10 au is a scale radius chosen to truncate the initial 
profile. The constant Eo is chosen such that the initial accretion rate at r ln is M = 5 x 10~ 8 M 0 yr -1 . 
With these choices, the initial disk mass is 0.05 M 0 , and the disk lifetime is around 8 Myr. These values 
are consistent with observations (Beckwith et al. 1990; Strom et al. 1989; Haisch, Lada & Lada 2001). 
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Solids in the disk are divided into two components, depending upon whether they are well-coupled to 
the gas (colloquially ‘dust’) or decoupled (‘planetesimals’). The dust component, of surface density Ed, is 
treated as a trace contaminant to the dominant gas component, with concentration C = Ed/E -C 1. This 
simplifies matters considerably, since the concentration obeys the equation (Clarke & Pringle 1988), 


dC 13/ BC\ dC 
~di ~ rZ V rS ) ~ Vr lh' 


( 5 ) 


where D is an effective diffusion coefficient and v r is the radial velocity. This equation allows for both 
advection of dust with the mean disk flow, and diffusion in regions of the disk where there is a gradient in 
the dust concentration. Although both are byproducts of turbulence, there is no reason why the diffusion 
coefficient should equal the viscosity 3 , but I assume that they have the same radial scaling and write D = £zq 
with £ a constant. Initially, I take C = 10~ 2 throughout the disk. The disk wind is assumed to entrain dust 
with it, leaving C in the material left behind unaltered. 

Coagulation of the dust will produce a population of larger (km sized and upward) bodies which are 
decoupled from the gas and follow Keplerian orbits. Eventually, these may form planets (e.g. Kenyon (2002), 
and references therein). It is generally believed that the rate limiting step in this process occurs during the 
early phases of dust settling (Weidenschilling 1980; Lissauer 1993), and that rapid growth occurs through 
the size range where radial drift is rapid. In these circumstances a local prescription for converting dust into 
planetesimals suffices. I adopt the simple fo 

dT, d 
dt 

dYip 

dt 

Here E p is the surface density in decoupled solid bodies, and tq is a characteristic time for dust settling at 
1 au (Nakagawa, Sekiya & Hayashi 1986). Observations suggest that significant dust evolution occurs on 
Myr timescales even at r > 10 2 au (Throop et al. 2001), which lends support to the theoretical conclusion 
that To is short. I assume a fiducial value To = 10 4 yr, but consider the effect of varying this parameter. 

If growing terrestrial planets attain masses of the order of M® prior to the dispersal of the gas disk, 
gravitational torques can once again couple the gas and solid components and drive rapid migration (Ward 
1997a). By using equation (6), which ignores this effect, I am therefore implicitly assuming that the final 
stages of terrestrial planet formation occur subsequent to the loss of the disk. In the terminology of Ward 
(1997b), this is the ‘slow’ mode of planet formation, which is often assumed to be appropriate for our own 
terrestrial planets (Wetherill 1990). The possibility of an alternate, ‘fast’ mode, in which Type I migration 
is an essential element, is briefly discussed in §3. 

Equations (1) and (5) are solved on a non-uniform radial grid, with 200 mesh points extending from 
0.067 au to 200 au. Standard numerical methods suffice for both equations, and can be validated using 
analytic solutions to special cases of both the surface density equation (Pringle 1981) and the contaminant 
equation (Clarke & Pringle 1988). Zero-torque boundary conditions are applied at the inner edge of the disk. 
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3 The relation between v and D could be determined via numerical simulations, though the properties of disk turbulence 
at the radii of interest remain subject to uncertainty due to the low ionization fraction and resulting suppression of the 
magnetorotational instability (Gammie 1996; Balbus & Hawley 2000). There are also constraints on the efficiency of mixing 
processes from the isotopic heterogeneity inferred for the Solar nebula (Lugmair Sz Shukolyukov 1998; Boss 2001). 
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3. Results 

Figure 1 shows the basic evolution of the aforementioned model. The initial disk has a gas surface 
density at 1 au that exceeds the minimum mass Solar nebula by a factor of about 5, and extends out to a 
few 10’s of au. With dust diffusion and planetesimal formation parameters £ = 1 and To = 10 4 yr, the build 
up of solids decoupled from the gas is rapid in the inner disk, but requires of the order of a Myr at 5-10 au. 
The final radial distribution of planetesimals is similar but not identical to the initial gas distribution. 

Once a massive planet has reached 1 au, the time remaining until merger is short. To a good approx¬ 
imation, therefore, the influence of giant planet migration on the distribution of solid material in the inner 
disk can be modeled by instantaneously sweeping the inner disk clear of planetesimals at some specified time 
^migrate- By varying fmigrate, we can study the extent to which dust continues to advect and diffuse to small 
radii, and thereby form new planetesimals. 

Figure 2 shows S p at r = 1 au as a function of ^migrate- The surface density of decoupled solid material is 
evaluated at the end of the calculation, after the gas in the disk has been lost. For the fiducial parameters of 
£ = 1, To = 10 4 yr, the final surface density of planet-forming material is heavily suppressed unless migration 
occurs extremely early - within the first 10 5 yr. Even taking into account the fact that, in the absence of 
migration, the final surface density of planetesimals is rather large (substantially in excess of the minimum 
mass Solar nebula value), it seems unlikely that a second generation of planetesimals would be formed in 
large enough numbers to permit terrestrial planet formation following migration. 

How sensitive is this conclusion to the rather uncertain input parameters of the model? Although the 
adopted t 0 = 10 4 yr has some theoretical support (Lissauer 1993), taking £ = 1 is simply a plausible guess. 
Since the concentration decreases toward small radii, where planetesimal formation is more rapid, more 
dust is able to repopulate the inner disk following migration if the efficiency of dust diffusion is enhanced. 
Increasing the timescale on which dust forms planetesimals obviously has the same effect. Figure 2 shows that 
by adopting £ = 5, or modestly increasing the planetesimal formation time to To = 3 x 10 4 yr, the timescale 
on which migration could occur without detrimental effects on subsequent terrestrial planet formation can 
be lengthened to perhaps a Myr. Since the disk lifetime is almost an order of magnitude in excess of this 
(Strom et al. 1989; Haisch, Lada & Lada 2001), migration occurring at a random epoch would in most cases 
still inhibit terrestrial planet formation using these values. This would not be true if To was as large as 
10 5 yr. However, with this To a single e-folding time for dust destruction at 5 au is in excess of a Myr. Early 
giant planet migration would then be unlikely, simply because there would probably not have been time to 
form a massive planet at 5 au via core accretion (e.g. Pollack et al. 1996) in the first place. 

Could low mass planets form further out in the disk, and subsequently migrate inward via gas interactions 
to fill the void at smaller radii? This is conceivable. The rate of Type I migration (in which the embedded 
planet fails to open a gap) is not limited to the inflow velocity of the gas, and can be extremely rapid for 
Earth mass planets at radii of the order of 1 au (Ward 1997a). If growing planets at somewhat larger radii 
attain such masses while the gas remains present, it is likely that they will suffer substantial radial migration, 
and pass through the terrestrial planet region. What is less clear is whether migration can be halted to leave 
terrestrial planets orbiting at small radii, given that the gas which drives migration will rapidly refill the 
inner disk following the passage of a massive planet. At very small radii, interaction with the star or its 
magnetosphere (Lin, Bodenheimer & Richardson 1996) provides a mechanism which can both halt migration, 
and allow further growth as additional low mass objects migrate inward (Ward 1997b). Further out, it is 
possible that the existence of sharp edges in the disk perhaps associated with the onset of efficient MHD 
turbulence within r ~ 0.1 au (Gannnie 1996) - could play a similar role (Kuchner & Lecar 2002). Additional 
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work is needed to assess the viability of these possibilities, and to determine how much gas would be accreted 
during the Type I phase. Initial numerical simulations by D’Angelo, Henning & Kley (2002) and Bate et 
al. (2002), which do not incorporate the complexities in the radial disk structure discussed above, find that 
both migration and accretion are very rapid for masses M ~ M®. This would suggest that Earth mass 
planets migrating through the disk would accrete significant gaseous envelopes. Simulations which include 
a treatment of the energy equation are needed to determine whether this is a robust result. 

Migration due to interactions with the gas disk is almost inevitable unless planets happen to form very 
close to the moment when the gas is dispersed. It is possible that further orbital evolution occurs subse¬ 
quently, driven by purely gravitational interactions with planetesimals (Murray et al. 1998) or other massive 
planets (Lin & Ida 1997; Ford, Havlickova & Rasio 2001; Marzari & Weidenschilling 2002). These processes 
could also impact the probability of terrestrial planet formation. In particular, if dynamical interactions are 
the cause of extrasolar planet eccentricities, then an unstable multiple planet system must be a common 
outcome of the planet formation process. The relevant question for terrestrial planet formation is then not 
which orbits are unstable in the present-day system, but which orbits were swept free of planetesimals during 
the early chaotic evolution of a massive planet system. This could be substantially more restrictive of the 
possibilities for habitable terrestrial planets than the constraints derived so far. 


4. Discussion 

The main technical result of this paper is that inward migration of massive planets, if it occurs more 
than ~ 10 6 yr after the formation of the protoplanetary disk, is likely to suppress the subsequent formation 
of terrestrial planets at radii of the order of 1 au. This suppression occurs because the gas that flows in to 
the inner disk, after the migrating planet has merged with the star, is depleted of dust as a result of already 
having formed planetesimals at larger radii. It is possible to choose parameters such that this conclusion can 
be evaded, but this requires either slow conversion of dust to planetesimals (in which case it is hard to see 
how the migrating massive planets could have formed, at least via core accretion), or very rapid diffusion of 
dust on a timescale an order of magnitude shorter than the disk’s viscous timescale. 

For future targeted searches of nearby stars for low mass planets, the direct implication is that stars 
with ‘hot Jupiters’ - massive planets at extremely small orbital radii of a ~ 0.1 au or less - may have been 
unable to form terrestrial planets, even though there are dynamically stable orbits for habitable planets 
in those systems (Jones, Sleep & Chambers 2001; Menou & Tabachnik 2002). Whether there are broader 
implications for the likely frequency of terrestrial planets depends first upon how common migration is, and 
second on whether there are stars that form terrestrial planets without forming any giant planets. The latter 
class of systems could certainly be common, though as their numbers are entirely unconstrained by existing 
observations any estimate is necessarily speculative. 

How common migration is among systems that form massive planets depends primarily on how fast 
migration occurs. If inward migration from a few au to merger with the star is rapid (compared to the disk 
lifetime), then it requires fortuitous timing to form a giant planet at a moment, close to the disk dispersal 
time, when it will not migrate significantly. Numerical simulations (Bate et al. 2002) indicate that the 
migration timescale (at least in the simplest disk models) is short for planets formed close to the snow 
line, and this is reflected in migration being the rule rather than the exception in statistical models of the 
extrasolar planet distribution (Armitage et al. 2002; Trilling, Lunine & Benz 2002). Among stars that either 
now possess massive planets, or once formed massive planets that were subsequently lost, migration through 
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to merger is likely to have occurred in the majority of systems (Trilling, Lunine & Benz 2002). As a result, 
it is possible that a high efficiency of massive planet formation could depress the frequency of terrestrial 
planets to a level substantially below the massive planet frequency. 


I thank Nick Schneider for posing difficult questions that prompted this work, and the referee for 
extremely helpful suggestions. 
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Fig. 1.— Build up of the surface density in planetesimals for a model with To = 10 4 yr and £ = 1. The 
dashed curves show results at log(t/yr) = 4, 4.5, 5, 5.5, 6 and 6.5. The upper solid curve shows the initial 
gas surface density distribution. 
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Fig. 2.— The final surface density in planetesimals at 1 au, as a function of when a single massive planet 
migrated through the inner disk. The solid curve shows the results for a model with to = 10 4 yr and £ = 1. 
The other curves each vary one parameter. The two dashed curves have longer planetesimal formation times 
(to = 3 x 10 4 yr, to = 10 5 yr), and the clot-dashed curve more efficient diffusion (£ = 5). The dotted line 
shows an estimate of the surface density in planetesimals at 1 au for a minimum mass Solar nebula model. 





